Numerical and Theoretical Study of a Monodisperse Hard-Sphere Glass Former 
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There exists a variety of theories of the glass transition and many more numerical models. But 
because the models need built-in complexity to prevent crystallization, comparisons with theory can 
be difficult. We study the dynamics of a deeply supersaturated monodisperse four-dimensional (4D) 
hard-sphere fluid, which has no such complexity, but whose strong intrinsic geometrical frustration 
inhibits crystallization, even when deeply supersaturated. As an application, we compare its behav- 
ior to the mode-coupling theory (MCT) of glass formation. We find MCT to describe this system 
better than any other structural glass formers in lower dimensions. The reduction in dynamical 
heterogeneity in 4D suggested by a milder violation of the Stokes-Einstein relation could explain 
the agreement. These results are consistent with a mean-field scenario of the glass transition. 

PACS numbers: 64.70. qj, 61.43.Fs,64.70.pm, 66.30.hh 



Though the conversation started generations ago, sci- 
entists still debate the nature of the glass transition [l| . 
The multiplicity of competing frameworks even led some 
to declare that glass theories are more numerous than 
glass theorists [2j. To be fair, a great deal of progress 
toward a microscopic description has been made over the 
last couple of decades, but many fundamental questions 
remain unanswered. The intrinsic complexity of many 
glass forming systems hinders progress. To avoid inter- 
ference from precipitous crystallization, simulated fluids 
require built-in complexity: bidisperse 0, 0], polydis- 
perse [f| or multicomponent Q mixtures, anisotropic Q 
or frustrated [8j pair interactions, etc. The situation is 
even more intricate in experiments, where good glass for- 
mers are, with only one reported exception 0], molecular 
liquids, polymer melts, or metallic alloys. Quantitative 
comparisons between the microscopic theories and these 
systems are formidably challenging, which impedes as- 
sessing the validity and limitations of the predictions [lfj| . 
A plain glass former, i.e., a monatomic, one-component 
liquid with an interaction potential as simple as possible 
would thus be greatly beneficial for the field's progress. 

In this Letter, we report numerical simulations of a 
glass former made of monodisperse 4D hard spheres (HS). 
Monodisperse HS have a single non-trivial thermody- 
namic parameter, the volume fraction r], which makes 
them the simplest model of isotropic fluids and crys- 
tals. Surprisingly, unlike in 2D and 3D, monodisperse HS 
in 4D and higher dimensions are stro ngly geometrically 
frustrated with respect to the crystal The simplex- 
based liquid structure is geometrically distinct from the 
crystalline order, which inhibits nucleation and facilitates 
glass formation 12, 13| • Compression studies of 4D HS 
suggest that the structural relaxation timescale becomes 
longer than the slowest accessible quenching rate around 
■q « 0.40 — 0.41 [l2|; 14 1 . A rough estimate gives the 



nucleation rate of reasonable system sizes for simula- 
tions to be eight to ten orders of magnitude slower than 
the structural relaxation timescale 13, 15 1. In monodis- 
perse 3D HS, in contrast, at rj « 0.54 multiple crystal 
nuclei form simultaneously on the structural relaxation 
timescale, while the onset of slow dynamics is generally 
agreed to be rj ~ 0.58. Bidisperse or polydisperse HS 
are thus used to study glass formation (e.g. [J, |5|). We 
perform molecular dynamics simulations with the event- 
driven package of Ref. [l2j in a system of N = 4096 par- 
ticles, in order to examine the 4D system's glass forming 
properties [lq . We expect finite-size effects to be small 
in this regime, as is the case at similar supersaturations 
in 3D for N = 512 0. The structure factor S{k) of 
the deeply supersaturated fluid remains liquid like at all 
densities, but the complete absence of crystallinity is also 
checked by an order parameter developed to detect nu- 



cleation [13j. Dynamically, the two-step growth of the 
mean-square displacement (MSD) (|rj(i)— r;(0)| 2 ), where 
Yi(t) is the position of the i-th particle, shows a lengthen- 
ing caging plateau with density, a signature of structural 
glass formers (Fig.[TJ). 

As an application of this simple model glass former, 
we consider the role of dimensionality in glass forma- 
tion, which is a subject of considerable theoretical dis- 
cussion [3, EM EH- The abrupt dynamical slowdown 
near the glass transition suggests the presence of a ki- 
netic and/or a hidden thermodynamic singularity. The 
most direct evidence for such a singularity is the growth 
of spatio-temporal fluctuations on the structural relax- 
ation timescale, which results in fast and slow moving 
regions in supercooled liquids. This dynamical hetero- 
geneity is typically monitored through four-point correla- 
tion function s [2011 and Stokes-Einstein (SE) relationship 
violations [2, l2l| . As for conventional continuous phase 
transitions, the impact of these fluctuations should be re- 
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duced in higher dimensions as the system becomes more 
mean-field-like. Comparing systems of different dimen- 
sionality should allow to better understand the glass tran- 
sition as a critical phenomenon and to test this mean-field 
scenario. A first attempt in this direction was recently 
made by Eaves and Reichman for a 4D binary Lennard- 
Jones (BLJ) model system (22[, but the complication of 
identifying the dimensional correspondence between BLJ 
systems makes quantitative comparisons difficult. 

Insights into the mean-field scenario of the glass tran- 
sition are gained by a MCT analysis of the dynamical 
results. A mean-field description of the thermodynamic 
"ideal" glass transition based on the replica theory sug- 
gests that HS undergo a dynamical transition at rj c , be- 
fore reaching the thermodynamic glass transition of the 
one-step replica symmetry breaking at t\k > r\ c 
This decoupling between dynamical and thermodynamic 
anomalies is demonstrated by simulation |24j . MCT is 
conjectured to be the dynamical counterpart of this ap- 
proach below rj c , because its mathematical structure is 
equivalent to the dynamical equations of a mean-field p- 
spin glass model 25[ for which the relation between the 
dynamical and the thermodynamic glass transitions is 
rigorously established [2(|. Though contentious, MCT 
is one of the most successful theories of the glass transi- 
tion. It uses static structural information, such as the ra- 
dial distribution function to provide first-principles 
predictions of the slow dynamics of fluids before their 
dynamical arrest [27|. At mild supercooling it qualita- 



tively captures the onset of the two-step decay of time 
correlation functions and the algebraic relaxation of the 
intermediate time regime. On approaching rj c it predicts 
a power-law divergence of the structural relaxation time 
instead of the well-known Vogel-Fulcher-Tammann be- 
havior This power law describes simulation and ex- 
perimental observations fairly well for a range of den- 
sities below the r\ c fitted from the dynamical data [27j |. 
but the divergence at rj c is rounded off due to activated 
events, which the theory does not capture. A failing of 
MCT is that it foretells a nonergodic freezing of the dy- 
namics at a much lower r\ c (or higher temperature T c for 
thermal systems) than the experimental and simulation 
glass transition point r/g (or T g ) Moreover the non- 
ergodic freezing point obtained by fitting the simulation 
data with the MCT power law are systematically lower 
than what the theoretical predicts 0, |2^. MCT also 
lacks an explanation for the violation of the SE relation, 
which is mostly attributed to strong dynamical hetero- 
geneity near the glass transition and is missing in the 
theory [3(J. Yet if MCT is indeed a dynamical mean- 
field theory, a reduction in the degree of heterogeneous 
dynamics by increasing dimensionality should improve 
the its agreement with simulation results. 

MCT is expressed as a series of nonlinear integro- 
differential equations for correlation functions such 
as the intermediate scattering function F(k, t) = 




FIG. 1: (Color online) (a) MSD and (b) F s (k,t) at 7] = 0.370, 
0.386, 0.389, 0.392, 0.395, 0.398, 0.401, 0.404, and 0.407, from 
left to right. Solid lines are MCT fits using e as control pa- 
rameter (see text). At the highest density r a deviates from 
the power law scaling of Fig. [2] 



N~ 1 {5p\ 5 _(t)Sp^(0)}, where Spk(t) is the density fluctu- 
ation in reciprocal space. Generalization of the theory to 
4D straightforwardly gives [3l[ 

F(k,t)+CllF(k,t)+ [ds M{k,t-s)F{k,s) = 0, (1) 



where Stf. = k B Tk 2 /mS(k), S(k) = F(k,t = 0) is the 
static structure factor, and M(k, t) is the memory kernel. 
M(k,t) can be further decomposed into fast and slow 
components M[k,t) = Mf ast (k,t) + Mmct(&,£) with 
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where V k {q,p) 
/167rfc 4 , c(Jfe) : 



q+k\ 

dpV k (q,p)F(q,t)F(p,t), 

\q-k\ 

(2) 

[f+c(q) + f-c(p)} 2 



7] ^4k 2 p 2 -fl 

{1 — l/S(k)}/p is the direct correla- 
tion function, p — N/V is the number density, and 
f± = k 2 ± (q 2 — p 2 ). Binary collisions dominate the fast 
decaying part of the memory kernel Mf ast (fc,i), which is 
conventionally determined by fitting simulation data [28j j . 
Here, the analysis is done for the self part of the interme- 



Jk.[vi{t)-Ti{0)]\ for 
which the MCT expression is similar to Eq. [2] 31]. Bc- 



diate scattering function F s (k,t) = (e l 



cause the MCT analysis is very sensitive to the details 
of the S(k) input, particular care is taken to interpolate 
and extrapolate (to larger k) the simulation data to a con- 
tinuous functional form. We extend the simulated g(r) 
for distances larger than half the simulation box with a 
damped oscillatory function. The Fourier transform of 
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FIG. 2: Power-law fit of D (A) and r^ 1 for k = 8.3 (•). Inset: 
D 1 / 1 and T a for k — 8.3, 7.5, and 4.2, from top to bottom, 
scale linearly with a shared intersection point. 



g(r) is in good agreement with the direct computation 
of S(k). We compute (see Fig. [T]) F s (k,t) for several 
wavevectors, one of which (k — 8.3) is close to the first 
peak of the structure factor and thus not too far from 
the dominant microscopic caging and relaxation length 
scales 1.32] . The decorrelation of density fluctuations dis- 
plays the characteristic caging plateau and allows for the 
extraction of the structural relaxation time r Q , defined 
as F s (k,T a ) = 1/e. 

MCT predicts (Eq. [J) that both r a and the diffusion 
coefficient D scale with the same power-law exponent 7 
and critical density rj c , i.e., r~ x ,D oc \r\ c — ?7| 7 , with 
7 MCT = 2.8 and r]f GT = 0.379. In Fig. [21 we fit the 
simulation data to this form by fixing 'y MCT and leav- 
ing T] c free. All the data except for the densest system 
follow a power law with rf^ = 0.407 and rf™ a = 0.406 
for all wavevectors. The discrepancy between 77^™ and 
ijl^l is smaller (0.2%) than for 3D HS (0.5%) 0], and 
the concordance with ?7^ CT slightly improves. It is in- 
teresting to note that r]l im is very close to the extrapo- 
lated dynamical arrest point inde pen dently obtained at 
very slow compression rates 



For consistency 
= 0.406, which 
all wavevectors. 



check, we also fit 7 while fixing rjl im 
gives 7 B m = 2.4 and 7?™ = 2.7 for 
In contrast, for 3D BLJ 0, [33|, which is the best char- 
acterized system, 7)5"™ ss 1.8 and 7^™ m 2.3, whereas 
7 MCT = 2.46. Another notable feature is that the power 
law fits over almost four time decades, reaching densities 
where e = l-r)/rfc im is less than 0.5% (Fig.©. In 2D and 
3D systems the structural relaxation timescale follow a 
power law for only two or at most three decades 0,0, 33- 
|, before the theoretical description breaks down be- 
cause of activated processes. 

We also consider how well MCT describes the full MSD 
and F s (k,t) curves, using e as input parameter to scale 
out the r] c dependence (Fig. []}. To tease out the long- 
time MCT behavior from Eq. [TJ the short-time F s (k,t) 
decay (up to t « 10°) is imported from the simula- 
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FIG. 3: (Color online) Violation of the SE relation in 4D 
monodisperse and 3D bidisperse ,4] and polydisperse [5J HS, 
where Dr a is normalized by its low density (-Dr Q ) re f value. 
MCT predicts only a very small violation of the SE relation 
for both 3D and 4D systems. 



tion curves. This procedure is equivalent to, but more 
direct than the standard fit of Mf ast (k,t) from simula- 
tion [28(. The concordance between simulation and cal- 
culated F s (k,t) is almost perfect for k — 8.3 and 11.9, 



up to densities where r a deviates from the power-law be- 
havior and MCT foretells the nonergodic transition. The 
calculated MSD also matches the simulated curves quite 
well, except for 77 = 0.404 and above. The correspon- 
dence between simulation and MCT is noticeably better 
than for lower-dimensional systems, where for every re- 
laxation curve both r\ c and k must be rescaled, in order 
to achieve a reasonable collapse 36]. Besides the power- 



law scaling, MCT successfully describes several other fea- 
tures in 3D, such as the time-temperature superposition, 
the von-Schweidler law at the beta relaxation regimes, 
the fc-dependence of the plateau height (the non-ergodic 
parameter), the exponents of stretched exponential re- 
laxation at the alpha relaxation regime j3(k), and the 



beta relaxation exponent b 27|. The concordance be- 
tween theory and simulation for all of these features is at 
least as good for 4D HS as for 3D systems. 

The quantitative improvement of MCT predictions in 
4D over 3D suggests that fluctuations are less important 
in higher dimensions. We analyze the system's dynam- 
ical heterogeneity through the SE relation to verify this 
hypothesis. Figure [3] shows Dr a for the full range of 
densities explored along with the corresponding 3D HS 
results 0, H| . In order to put the different dimensions on 
an equal footing, we plot Dr a against the diffusion coef- 
ficient D. The SE relation holds for diffusivities over a 
decade smaller in 4D than in 3D (up to e ~ 1%), which is 
slightly before D and r a start deviating from the MCT- 
predicted power law in Fig. [2] If a similar degree of poly- 
dispersity were used in 4D as in 3D we expect the SE 
violation to be suppressed even more [37]. The suppres- 
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sion of the SE relation violation and the better agreement 
of all available dynamical observables with MCT predic- 
tions suggest that 4D HS are dynamically more mean- 
field-like than the 2D and 3D equivalents. The improve- 
ment of the agreement with MCT is however incremental, 
which hints that if an upper critical dimension d c for the 
glass transition exists, it is larger than four. This result 
is consistent with general theoretical arguments that give 
d r = 8 18, 381. Biroli et al. also obtained that the SE 



relation violation should scale as Dr a ~ e 



d/4-2 



IB. The 



growth in Fig. [3] is not inconsistent with this scaling, but 
the SE relation violation for the range of densities we ex- 
plore is too mild to be conclusive about the exponent. A 
more direct measure of dynamical heterogeneity would be 
to compute the four-point correlation function Xi{t)i but 
the lack of statistical accuracy and the relatively small 
system sizes prevent us from reporting the results here. 
We will consider Xiif) lTl future studies. 

The 4D monodisperse HS fluid we study is quite con- 
venient to examine the glass transition. Its simplicity 
and slow nucleation rate allow high-accuracy compar- 
isons of its glass- forming properties with microscopic the- 
ories. The agreement of the system with MCT, which is 
broader than for any lower-dimensional equivalents, and 
the strong suppression of SE relation violation are consis- 
tent with the dynamical mean-field scenario of the glass 
transition [l4j| . The results also suggest that 4D is still 
below the upper critical dimension, if it exists, because 
the fluctuations due to activated processes round off the 
sharp dynamical singularity. Given that even the mean- 
field picture has not been rigorously established for the 
structural glass transition, the study of clean, higher- 
dimensional systems such as 4D HS is likely to play a 
crucial role in assessing the validity and limitations of 
the various glass theories. 
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